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O ' Abstract 
(N 

Low-Q 2 photons do not resolve partons in the proton, which gives problems when 
applying the deep inelastic scattering formalism, such as an unphysical, negative 
gluon density extracted from data. Considering instead hadronic fluctuations of the 
| photon, we show that the generalised vector meson dominance model (GVDM) gives a 

good description of the measured cross section at low Q 2 , i.e. reproduces F 2 (x, Q 2 ), 
using only few parameters with essentially known values. Combining GVDM and 
parton density functions makes it possible to obtain a good description of F 2 data 
over the whole range of x and Q 2 . 

<N 
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O ■ 1 Introduction 

Experimental measurements on electron-proton (ep, and also fip) scattering are usually 
Ph! interpreted in terms of the theoretical formalism for deep inelastic scattering (DIS). The 
differential cross section is then expressed in terms of proton structure functions given by 
the density functions for different partons, i.e. q(x, Q 2 ) and g(x, Q 2 ) for quarks and gluons 
carrying a fraction x of the proton's energy-momentum when probed with the scale Q 2 . 
The structure function F 2 , which gives the dominant contribution to the cross section, is in 
leading order given by F 2 (x,Q 2 ) = ^2 q £q {xq(x,Q 2 ) + xq(x, Q 2 )) while the gluon density 
enters indirectly via the logarithmic Q 2 dependence of perturbative QCD. 

This formalism has also been applied to F 2 data at low photon virtuality Q 2 , where the 
exchanged photon is not far from being on-shell. Parametrising such F 2 data in terms of 
quark and gluon density functions results in gluon distributions that tend to be negative 
for small x at small Q 2 (e.g. x ~ 10~ 4 , Q 2 ~ 2GeV 2 ) E]. The reason for this is 
that the DGLAP evolution, driven primarily by the gluon at small x, otherwise gives too 
large parton densities and thereby a poor fit to F 2 in the genuine DIS region at large 
Q 2 . Although one may argue that the gluon density is not a directly observable quantity 
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and hence might be negative, it certainly is in conflict with the interpretation of the 
probability for a gluon with momentum fraction x in the proton. In particular, such 
a gluon distribution could be just an effective description for a more proper theoretical 
understanding. It need not have the same universality as proper parton density functions, 
thus giving incorrect results when applied to other interactions. For example, differences 
in the predicted Higgs production cross section (dominated by gg — ► H) at the Tevatron 
and LHC arise depending on whether the gluon parametrisation is forced to be positive 
definite or allowed to be negative at small x [2]. 

In this Letter, we argue that the root of the problem is the application of the formalism 
for DIS also in the low-Q 2 region, where the momentum transfer is not large enough that 
the parton structure of the proton is clearly resolved. The smallest distance that can 
be resolved is basically given by the momentum transfer of the exchanged photon through 
d = 0.2/ y^Q 2 , where d is in Fermi if Q 2 is in GeV 2 . This indicates that partons are resolved 
only for Q 2 > 1 GeV 2 . For Q 2 < 1 GeV 2 , there is no hard scale involved and a parton basis 
for the description is not justified. Instead, the interaction is here of a soft kind between 
the nearly on-shell photon and the proton. The cross section is then dominated by the 
process where the photon fluctuates into a virtual vector meson state which then interacts 
with the proton in a strong interaction. This is the essence of the vector meson dominance 
model (VDM), for a review see 

In the following we use the original generalised vector meson dominance model (GVDM) 
[3] for ep scattering at low Q 2 . We show that it gives a good description of the recent 
HERA data extending the Q 2 region to very low values, which are of particular impor- 
tance for the GVDM approach (for a review of GVDM models, see [S]). Furthermore, the 
GVDM model based on hadronic fluctuations of the photon is natural to combine with 
our model [6 j for hadronic fluctuations of the target proton, which has been used to derive 
the non-perturbative x-shape of the proton's parton density functions. Combining parton 
density functions including DGLAP evolution 7J with GVDM gives a good description of 
data over the full Q 2 region. This extends earlier work [HJ E] on applying GVDM and is 
complementary to theoretical developments where GVDM is connected with a QCD dipole 
approach 011111121113. 

2 Vector meson dominance model for ep at low Q 2 

The occurence of quantum fluctuations implies that a photon may also appear as a vector 
meson such that the quantum state should be expressed as 

l7) = C |7o) + ^f \V)+ ! dm(...) (1) 

y JV J mo 

The first vector meson dominance model included only the sum over the vector meson 
states V = p°,u,</>. . ., whereas the generalised model |lj also includes the integral over a 
continuous mass spectrum (not written out explicitly in eq. (Q)). 
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This hadronic fluctuation of the photon then interacts with the target proton with a normal 
hadronic cross section dominated by soft processes without any hard scale involved. Total 
cross sections for different beam hadrons at different energies are well measured and given 
by standard parametrisations to be discussed below. The overall cross section is then a 
convolution of the photon-to-meson fluctuation probability with the meson propagator and 
the meson-proton cross section. 

In ep scattering 3 data is given in terms of the proton structure function F2 extracted 
from the differential cross section da/dxdQ 2 for electromagnetic interactions (one-photon 
exchange), since the weak interactions are completely negligible for Q 2 -C rn 2 zw . The 
structure function F 2 can be expressed as [TH |3] 

in terms of the total cross sections and for transverse and longitudinal virtual 
photons. 

These cross sections are obtained by squaring the amplitude involving expression (0) whose 
continuous part results in a double mass integral 

j m 2 ^ m2 ^ m/2 ^ T '(m^(y)('m /a +Q 2 7 Off-diagonal contributions having m 7^ m! [12] are 

normally neglected in phenomenological studies on nucleons, although they cannot be ne- 
glected for nuclei [HH ITo]. Since we here only consider nucleons, we take this integral to be 
diagonal, i.e. pT 1 L(W 2 ,m 2 J m' 2 ) = p TjL (W 2 ,m 2 )5(m 2 — m' 2 ). The spectral weight function 
Pt(W 2 ,m 2 ) is phenomenologically chosen to fit data, e.g. pt = m 2 ] /m A to obtain scaling 
at larger Q 2 , while pL = ic^zPT- In this GVDM approach, the resulting cross-sections are 
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In the sums over the discrete vector meson states one recognises the well-known factors 
Ana/ fy (involving the vector meson decay constant fy) which give the probabilities of the 
fluctuations 7 — > V for real photons, followed by the squared propagator of the meson 
with mass my and the meson-proton total cross section oy v . The terms proportional 
to crop = r c o~-yp (defined exactly below) originates from the integral over the continuous 

3 The DIS variables are defined through Q 2 = -q 2 = —{p e — p' e ) 2 , x — Q 2 /2P ■ q, y = P ■ q/ P ■ p e in 
terms of the four- momenta P,p e ,p' e ,q of the incoming proton, incoming and scattered electron and the 
exchanged photon, respectively. 
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vector meson mass spectrum with a lower limit given by the parameter m . The parameters 
4V — (J vp/ 'vp an d £c — a cpl a Cp accounts for the possibility of different cross sections for 
transverse and longitudinal polarisation states. It is assumed that they are independent of 
x and Q 2 and expected that they are less than unity. 

The total cross-sections oy p and be directly taken as the well known and generally 

used parametrisation [T^J 

a{ip -> X) = A iS e + (5) 

for the total cross section of a particle i on a proton. The first term is for pomeron 
exchange and the second one for reggeon exchange. The energy dependence is given by the 
parameters e « 0.08 and rj « 0.45 which are universal and obtained from fits to a wealth 
of data on total cross sections, whereas the normalisation parameters Ai,Bi are different 
for different particles. At high energies the reggeon term can be neglected in comparison 
to the dominating pomeron term. 

This parametrisation applies not only to the vector mesons (i = V) but also to photons 
(i = 7) which are on-shell or nearly so. Thus we have oy v = Ays e + Bys~ v and o lv = 
A 1 sl i + B^s' 71 . The fractions of the cross section accounted for by the discrete vector 
mesons V are then ry = y^- ■ and we can specify tq = 1 — XV r v as the fraction from 
the continuous mass spectrum. 

Inserting these GVDM expressions for ot,l i n ec i- © one obtains 

^ (l-x)Q 2 f f m 2 v V f ^ Q 2 



+ r c 



t s HI, 

4c 



Q 2 



^ + ^ln(l + ^)lU 7 ^ (6) 



where the following approximations, which are justified for the region of x and Q 2 of HERA 
data, have been made: In the prefactor the term 4x 2 m 2 /Q 2 <C 1 and is hence neglected. 
The last factor originating from <jy p and Oc v only includes the pomeron term, since the 
reggeon term is negligible, and the energy variable is s 1P = Q 2 ^~^+ m p ~ Q 2 /x at small-x. 

The parameters involved in eq. © are all essentially known from GVDM phenomenology. 
The values ry =p ^^c — 0.67,0.062,0.059,0.21 are quite well determined 0. Although 
m ~ 1 GeV is expected P , it is not well known and is here taken as a free parameter. The 
parameters £y are assumed to be the same for V = p, and expected to be £y ~ 0.25 
based on the early study in [I] and supported by [13] including recent HERA data. A 
similar magnitude is expected for £c- Lacking established numbers and wanting to have as 
few parameters as possible, we use the common parameter £ = £y — £c as a free parameter 
to be fitted. For the pomeron intercept parameter the value e = 0.09 has been obtained in 
recent fits ^H|, but we take it as a free parameter in order to check the expected consistency 
with this universal value. Also the overall normalisation constant A 1 of the photon-proton 
cross section is taken as a free parameter. Thus, we have the four parameters £, m , e, A^ 
to be fitted to data. 
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3 Comparison to F 2 data 
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Figure 1: F 2 at low Q 2 : HERA ep data from ZEUS [201 compared to GVDM as in eq. © 
(full curves). Model results are also given when the longitudinal contribution of the con- 
tinuum is excluded (£c* = 0) and when excluding the continuous contribution altogether 
(setting rc = 0) giving VDM. 



The GVDM expression for F 2 in eq. flOJ) gives a very good description of the HERA data 
on i<2 at low Q 2 , as shown in Fig. ^ The fit gives x 2 /d.o.f. = 87/(70 — 4) = 1.3 with 
parameter values as expected: e = 0.091, £ = 0.34, m = 1.5 GeV just above the discrete 
vector meson masses and A 1 = 71 fib in accordance with the measured photon-proton cross 
section (cf. ^Hj)- This demonstrates that for Q 2 clearly below 1 GeV 2 the HERA ep cross 
section can be fully accounted for by GVDM using parameter values as determined from 
old investigations related to fixed target data. 

For completeness, both the transverse and longitudinal contributions to the integral over 
the continuous mass spectrum are here included, although the latter is numerically small 
as demonstrated in Fig. ^ VDM, which lacks the continuum part, falls below the data 
and decreases too fast with Q 2 . This Q 2 behaviour becomes even worse if the longitudinal 
contribution is neglected (i.e. £y = 0), as is done in some simplified treatments of VDM. 
The Q 2 dependence of these different contributions is shown in Fig. |2j 
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Figure 2: The Q 2 dependence of F 2 from GVDM (full curve) with its contributions from 
transverse (T) and longitudinal (L) parts of the discrete vector meson spectrum (VDM) and 
the continuous (Cont.) mass spectrum. Data from SLAC [21] are included for comparison. 



We have also compared with data on F 2 from SLAC [21] and NMC [221- Due to the 
lower energies of these fixed target experiments, one must here include also the reggeon 
term in the Donnachie-Landshoff parameterisation of the total cross section and we use 
r] = 0.45, £> 7 = 90 fib (cf. [inillll)- Keeping the values of the other parameters fixed, we 
obtain good agreement as long as x and Q 2 are not too large (cf. j2Hj)- 

At larger Q 2 , this original GVDM does not have the correct behaviour since F 2 in eq. (jB]) 
increases with Q 2 for all x. This can be cured phenomenologically by introducing for the 
spectral weight function mentioned above a suitable form = N In {W 2 /am 2 )/m i [H]. 
With suitable values of the free parameters m , N, a it is then possible to reproduce HERA 
F 2 data also at larger Q 2 . A theoretically more advanced alternative is to instead include 
off-diagonal contributions ^Hl ^] . This connects naturally to the dipole formalism of DIS 
and include effects of perturbative QCD evolution. This off-diagonal GVDM framework 
should then apply in the full Q 2 region, as long as x is sufficiently small, and HERA data 
can here be reproduced [T2"j . 

At high Q 2 the conventional description is in terms of parton density functions, which 
also includes the large- a; valence region. As argued above, this approach does not apply 
at very small Q 2 and one must therefore complement it with GVDM to account for this 
region. To cover the full x and Q 2 region one should combine these two descriptions, but 
due to the confinement problem, there is no proper theoretical way to do the transition 
from GVDM formulated in a hadron basis to the parton model in a parton basis. Although 
GVDM can be extended to large Q 2 , this would imply double counting if combined with the 
conventional parton description. To use the latter one must, therefore, phase out GVDM. 

Thinking in terms of the resolution scale discussed above, it is quite natural that the 
original hadron-based GVDM only applies at low Q 2 and there should be a transition to 
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the DIS formalism of resolved partons at high Q 2 . In particular, the total cross sections 
°Vp5 a Cp used in GVDM applies to soft hadronic processes for (nearly) on-shell particles. 
It is therefore very reasonable to phase out GVDM at larger Q 2 by applying a form fac- 
tor suppression. A factor like m v /(m v + Q 2 ) [21| would, however, ruin the very good 
description at low Q 2 seen in Fig. ^ Instead, a sharper transition to DIS in the region 
Q 2 = 0.6 — 1.5 GeV 2 is required. This is in accordance with the rather abrupt change of the 
slope parameter A in F 2 (x) ~ x~ x observed in HERA data at Q 2 pa 1 GeV 2 [23] and may 
be seen more generally as a rather sharp transition from soft, non-perturbative to hard, 
perturbative QCD dynamics. 

We therefore introduce the phenomenological form factor (Qc/Q 2 ) a for Q 2 > Q 2 C to phase 
out GVDM above a critical Q 2 C . As shown in Fig. El a good description of HERA F 2 
data at intermediate Q 2 can then be obtained by combining GVDM and parton density 
functions that fit HERA F 2 data at larger Q 2 . This requires Q 2 C ~ 1 GeV 2 as expected from 
the discussed transition, and a pa 2 giving ~ Q~ A as a reasonable form factor damping. 
The exact values of the parameters are fitted and depend on the details of the DIS parton 
densities. With such a form factor suppression, the GVDM contribution is negligible for 
Q 2 > 4 GeV 2 (see fig. El), where DIS parton density parametrisations are usually considered 
trustworthy. Any parametrisation of parton densities which is good enough to reproduce 
the measured F 2 in the DIS region can be used, provided the GVDM component is taken 
into account when low-Q 2 data are included in the fits. 
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Figure 3: F 2 at intermediate Q 2 : contribution of GVDM with a form factor (1.24/Q 2 ) 1 ' 63 
(full curve) and the complete model (dashed curve), including also DIS parton density 
functions from our model, compared to HI data 26J. 

For Fig. El we have, however, used a physically motivated model |Hj where the parton 
momentum distributions are obtained from gaussian fluctuations having widths related to 
the uncertainty relation and the proton size. Valence distribtuions arise from the 'bare' 
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proton, whereas sea distributions originate from mesons in hadronic fluctuations of the 
proton \p) = a \p ) + otp^Poit ) + a nn \mi + ) + . . . + aAK\AK + ) + . . . This gives the x-shape 
of the parton densities at ~ 1 GeV 2 and the DGLAP equations are then used to evolve 
to larger Q 2 , resulting in a good fit to HERA F% data using only six parameters with 
physically motivated values 0. Furthermore, this model gives |27 u v (x) 7^ d v (x) and 
u(x) 7^ d(x) in qualitative agreement with data, as well as s(x) 7^ s(x) of interest for the 
NuTeV anomaly [Sl- 
it is interesting that combining these models involving quantum fluctuations of both the 
photon and the target proton results in a good description of the ep cross section, or 
equivalently F 2 , at both low and high Q 2 . 

4 Conclusions 

The conventional parton model formulation of deep inelastic scattering is not applicable 
at very low Q 2 , where no hard scale is available to resolve the partons. Instead, HERA F 2 
data are here well reproduced by the original generalised vector meson dominance model, 
including contributions from a continuous mass spectrum and longitudinal polarisation 
states, and using parameter values in agreement with old analyses at fixed target energies. 
At large Q 2 , GVDM with off-diagonal contributions can be used as long as x is small. To 
cover the full x-region, including the valence part, the proton structure must be introduced 
via parton density functions in the conventional DIS formalism. 

We have shown that one can combine the GVDM and parton density descriptions in a 
two-component phenomenological model. GVDM then accounts fully for the cross section 
below Q 2 < 1 GeV, but although it contributes also at large Q 2 it must here be phased 
out in order to avoid double counting with the standard parton density formulation. We 
have found that a form factor damping of GVDM gives a smooth transition into the deep 
inelastic region described by parton distribution functions. Here, any good parametrisation 
of parton densities can be used, provided the GVDM component is taken into account at 
low Q 2 as shown above when fitting the parameters. In this way one obtains a good overall 
result at both low and high Q 2 . In particlar, there is no need for a negative gluon density 
in the region of low x and low Q 2 . The reason is that the cross section is here dominated by 
the GVDM contribution, which is based on fundamental quantum fluctuations that should 
not be neglected. 

Acknowledgments: We are grateful to Dieter Schildknecht for interesting and helpful 
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